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We investigate the two-dimensional Hubbard model with next-nearest-neighbor hopping, t' , using 
the dynamical cluster approximation. We confirm the existence of a first-order phase-separation 
transition terminating at a second-order critical point at filling nc{t') and temperature Tpsit'). We 
find that as t' approaches zero, Tps{t') vanishes and nc{t') approaches the filling associated with the 
quantum critical point separating the Fermi liquid from the pseudogap phase. We propose that the 
quantum critical point under the superconducting dome is the zero-temperature limit of the line of 
second-order critical points. 



INTRODUCTION 

Strongly correlated electronic materials, which include 
high-temperature superconductors, heavy fermions, and 
magnetic compounds, are characterized by competing 
phases and complicated phase diagrams. These compet- 
ing phases can lead to quantum criticality when one of 
the transition temperatures is driven to absolute zero as 
a function of a non-thermal control parameter such as 
pressure or dopingii^. While the physics of a conven- 
tional phase transition is driven by thermal fluctuations, 
near a quantum critical point (QCP), quantum fluctua- 
tions affect the properties of a material up to surprisingly 
high temperatures^. In particular, transport measure- 
ments of hole-doped cuprates suggest the presence of a 
QCB^ lying beneath the superconducting (SC) dome^ 
Although it is believed that this QCP dominates the 
phase diagram, its nature is still unknowi*^ with compet- 
ing scenarios emphasizing the role of bosonic or fcrmionic 
fluctuations 1^ 

In this work, we provide evidence for the nature of the 
QCP in the Hubbard model of the cuprates by a system- 
atic study of its phase diagram. Our results suggest that 
the QCP is not due to order in the pseudogap (PC) re- 
gion, but rather is the zero-temperature limit of a line of 
second-order critical points associated with a first-order 
phase-separation transition (see Fig. [T| . The control pa- 
rameter for this transition is the next-near-neighbor hop- 
ping parameter, t' . 

Although it is possible to have a QCP not associated 
with any obvious order parameter, as in the case of a 
localization transition, in most QCPs a continuous order 
parameter vanishes at T = for a particular value of the 
controlling energy scaled. Less common is a QCP asso- 
ciated with a first-order transition, but this is possible 
when the first-order transition terminates at a second- 
order critical point which is driven to zero by tuning the 
relevant parameter. For example, in Sr3Ru207-i^, as a 
function of the field angle, a first-order meta-magnetic 
transition is driven to T = yielding quantum critical 




FIG. 1. (Color online) Schematic phase diagram of the 2D 
Hubbard model in the temperature (T), chemical potential 
(/i), and next-near-neighbor hopping {t') space. For t' > the 
first-order phase-separation transition terminates at a second- 
order critical point at doping Uc and temperature Tpa. The 
line of second-order critical points {Tps,nc), approaches the 
QCP on the t' = plane. This is the critical point separating 
the pseudogap (PC) from the Fermi-liquid (FL) region. 



phenomena. 

The latter case is consistent with the scenario of Fig. [TJ 
where the role of the field angle is played by t' . For pos- 
itive t' and below a certain temperature, the system un- 
dergoes a first-order phase-separation transition. In this 
region, two solutions with diff'erent densities coexist for 
a given chemical potential, More recently, this 

first-order metal-insulator transition was studied by Gull 
et ali^ and used to map out the phase diagram of this 
model in the space of interaction strength and t'. Since 
these two phases have the same symmetry, this transition 
terminated in a second-order critical point at tempera- 
ture Tps and critical filling Uc where the charge suscep- 
tibility diverges. By increasing t' , T^g increases and the 
critical point becomes numerically accessible. For low t', 
it is no longer accessible, but its presence is evidenced by 
a peak in the charge susceptibility. 

Our starting point is the two-dimensional (2D) Hub- 
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bard Hamiltonian, 



■ U E ni-fTiii, 



(1) 



where = —2t (cos + cos ky)—4t' (cos fc^; cos ky — 1) is 
the tight-binding dispersion as a function of the hopping 
t between nearest neighbors and t' between next-nearest 
neighbors, cjj.^(cko-) is the creation (annihilation) opera- 
tor for electrons of wave vector k and spin a, rii^ = cl^Ci„ 
and U is the on-site Coulomb repulsion. 



METHODOLOGY 

We solve the Hubbard model within the dynamical 
cluster approximation (DCA)ii on a A^c-site cluster. The 
DCA is a cluster mean-field theory which maps the orig- 
inal lattice model onto a periodic cluster of size Nc = 
embedded in a self-consistent host. Spatial correlations 
up to a range Lc are treated explicitly, while those at 
longer length scales are described at the mean-field level. 
However the correlations in time, essential for quan- 
tum criticality, are treated explicitly for all cluster sizes. 
To solve the cluster problem, we use weak-coupling ex- 
pansion continuous time quantum Monte Carlo (QMC) 
methodi^ii^ with highly optimized blocked and delayed 
updatesii, a determinant QMC method which scales lin- 
early in the inverse temperature^^, as well as Hirsch-Fye 
QMGi^iSa. The fast determinant QMC was used to ob- 
tain a converged DCA solution, Hirsch-Fye QMC was 
used to calculate lattice susceptibilities, and continuous 
time QMC was used as a control for systematic error. 
The unit of energy is t in the entire paper. 

To make contact with previous results, we perform sim- 
ulations with [/ = 6, but we find that the phase separa- 
tion becomes more prevalent for U — & for which we 
present most of our results. We calculate the filling, n, 
versus /i and the compressibility, dn/d^, by taking the 
numerical derivative. We also calculate various suscepti- 
bilities including the charge Xc(Q = 0, T), spin and pair- 
ing susceptibilities by solving the lattice Bethe Salpeter 
equation using the renormalized DCA vertices^S. Note 
that as a consequence of the fluctuation-dissipation the- 
orem, Xc(Q — 0,T) is identical to dn/dfi, but we keep 
both terms, compressibility and susceptibility, to identify 
the method by which they are calculated. 



RESULTS 

In Fig. [2Ka), we plot n versus fj. ior U = 6, T = 0.077 
and different values oit' , ranging from 0.0 to 0.4, on a 16- 
site cluster. The filling increases monotonically with the 
chemical potential and shows a pronounced flat region, 
associated with the Mott gap, especially for t' < 0.4. 
The most interesting feature of n{fi) is an inflection ap- 
parent at flnite doping, which becomes more pronounced 




FIG. 2. (Color online) Filling, n (solid lines), and compress- 
ibility, dn/djj, (dashed lines), plotted vs chemical potential, fi, 
for various values oft' for (a) U = 6, Nc = 16, and T — 0.077, 
and (b) U = 8, Nc = S at different temperatures. The unit 
of energy is t in all figures. The critical filling, where the 
compressibility peaks, is plotted in the corresponding inset. 
In (a) when t' — >■ the peak in the charge susceptibility is 
located at the QCP identified previouslj". 



for larger values of t'. The inflection in n(^) translates 
into a peak in the compressibility (shown on the right 
axis). The peak becomes sharper and moves closer to 
half filling as t' increases. For t' > 0.3, the peak disap- 
pears, as does the plateau in n(fj.) near half filling, associ- 
ated with the gap. The value of the critical filling at the 
peak, ric, versus t' is plotted in the inset. Note that for 
t' = 0, ric = 0.86 agrees with the filling 0.85 of the QCP 
identified previously for these paramcters^i*^ and a crit- 
ical filling separating two Fermi-liquid (FL) regions in a 
closely related t — J modeli^ These results suggest that 
the QCP may be associated with charge fluctuations. 

To explore this association, we study the behavior of 
the bulk charge susceptibility, Xc(Q = 0,T) and its di- 
vergence as t' — > 0. Unfortunately, for this cluster and 
parameters, the minus sign probleni^i limits our ability 
to access temperatures low enough to sec a divergence 
in the charge susceptibility. The minus sign problem be- 
comes worse when the cluster size, U or increases or 
the temperature decreases. However, the cluster size, in- 
teraction U and t' affect the phase diagram in different 
ways. In previous studiesii, we found that for t' = 0.3, 
clusters with Nc — 8, 12 and 16 have roughly the same 
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FIG. 3. (Color online) Inverse bulk charge susceptibility vs 
temperature when U = 8, Nc — 8 for several values of t' . The 
values of the critical filling Uc shown in the legend correspond 
to the maximum of low temperature compressibility, or the 
filling where it first diverges. 



strengthens the Mott gap for this t' as can be seen in 
the persistence of the flat region in near n = 1 for 
t' = 0.4 [Fig.Hb)]. 

The charge fluctuations associated with phase- 
separation influence the SC phase diagram. This is shown 
in Fig. H for Nc = 8, U ^ 8, and t' = 0.0, 0.1, and 0.3. 
The pscudogap temperature, T*, obtained as the tem- 
perature where the bulk spin susceptibility peaks (see 
Ref. [2lh . is also plotted. For t' = 0, T* vanishes at 
the QCP; which for this smaller cluster has moved to 
ric = 0.88. Note that for t' = 0, the SC dome is centered 
on the QCP, suggesting that superconductivity is asso- 
ciated with the quantum fluctuations. For t' = 0.1 and 
t' = 0.3, the SC dome contains the point where T* — s- 0, 
but is not centered around Uc- Instead, the second-order 
point is found on the low-doping side of the SC dome. 
Note that the maximum SC transition temperature in- 
creases slightly with increasing t' which is in agreement 
with previous results for a four-site cluster, and the same 
interaction strengtb^^. 



phase-separation transition temperature when U ~ S. 
Here, using the same interaction strength, we find that 
the cluster size effects are stronger for smaller values of 
t'. As the cluster size is decreased, the charge suscep- 
tibility is somewhat suppressed, and more significantly, 
the critical doping moves towards half filling. With in- 
creasing U in the range from [/ = 4 to 8, the peak in the 
charge susceptibility moves to lower fillings and higher 
temperatures. So, despite the worse minus sign problem 
associated with larger values of the interaction, the re- 
gion of divergent charge fiuctuations becomes larger and 
more accessible for U = 8. For this reason, from this 
point on we will use a smaller cluster with Nc = 8 and 
U = 8, in order to access the second order critical points 
and investigate their relationship to superconductivity. 

In Fig. [2lb), wG plot n and dn/d^ versus /i for sev- 
eral t' and temperatures. Similar to Fig. [D^a), we see a 
cusp emerge in the compressibility. As the temperature 
is lowered, the peak in the compressibility is enhanced for 
t' = 0.1 and 0.2. However, for a larger t' , e.g., t' = 0.3, 
and for T < 0.1, wc find hysteresis between two stable 
solutions with different values of n for the same value of 
/iii. The presence of hysteresis indicates that the system 
has undergone a first-order phase-separation transition. 

We explore the line of second-order critical points of 
these first-order transitions as t' changes using the charge 
susceptibility as shown in Fig. [31 Here, the inverse charge 
susceptibility at Uc is plotted versus temperature for dif- 
ferent values of t' when [/ = 8 and Nc = 8. The critical 
filling identified in the legend is determined as the filling 
where the compressibility either diverges, or is peaked 
at the lowest accessible temperatures. We find that the 
temperature of the second-order critical point increases 
with increasing t' and that it moves towards half filling. 
However, unlike the J7 = 6 results shown in Fig. [2Ia), 
the critical point appears to avoid half filling even for 
t' = 0.4. The stronger Coulomb interaction U = 8 also 



DISCUSSION 

A detailed study of the phase diagram of the 2D Hub- 
bard model with next-ncar-neighbor hopping has allowed 
us to identify the nature of the QCP under the SC dome. 
We argue that QCP is the terminal point of a line of 
second-order critical points associated with first-order 
phase-separation transitions. The critical temperature 
is driven to zero as t' 0. For positive t' , a Mott liq- 
uid and a Mott gas coexist at fixed In real mate- 
rials other parameters, such as electron-phonon interac- 
tion, inter-site electron-electron interaction or inhomo- 
gencities, might play a role similar to t' . 

Generally, it is accepted that the model with t' > 
describes the electron-doped and that with t' < de- 
scribed the hole-doped cuprates. We find that the model 
for t' > Q does not display quantum criticality, but rather 
classical criticality. The QCP is found only for i' = and 
as is known from other quantum critical systems, it will 
strongly affect the system for a wide range of parameters 
and temperatures around this point, including t' < 0, the 
model for the hole-doped cuprates. 

The relationship of the superconductivity with the 
QCP at t' = is not yet clear, but the fact that the dome 
is centered at the QCP suggests that the incipient phase 
separation creates conditions favorable for superconduc- 
tivitj*2S. The Mott liquid phase may provide regions 
where the spin-mediated pairing interaction is strong and 
the Mott gas may provide regions where there are quasi- 
particles to pair. It may also be that the incipient charge 
fluctuations when combined with the antiferromagnetic 
spin fluctuations enhance the pairing within a narrow 
region near the QCP~. The phase-separated region 
might also be related to the pervasive inhomogeneities 
observed in cuprates which led to theoretical scenarios 
for an inhomogeneity-bascd pairing mechanis m^^'^^ or 
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FIG. 4. (Color online) SC transition temperature Tc (solid 
lines), pseudogap temperature T* (dotted lines), and critical 
points (stars) for the 2D Hubbard model with U = 8, Nc = 8, 
and t' = 0, 0.1, and 0.3. For all t' the T* line terminates inside 
the SC dome. For t' — 0, this termination point coincides with 
the QCP of the phase-separation transition (star in the left 
panel). For t' = 0.1 and 0.3 the second-order critical point 
is at critical filling ric = 0.93 and 0.97, respectively, which is 
above the SC optimal filling. 



CONCLUSION 

We find that the QCP at i' = of the 2D Hubbard 
model is the zero-temperature limit of a line of second- 
order critical points associated with a first-order phase 
separation transition, which occur at finite temperature 
when t' > 0. The filling associated with the second-order 
critical point is determined from the peak position in 
the compressibility versus chemical potential. The peak 
grows by decreasing the temperature or increasing the 
cluster size or the interaction strength. We also show 
that for t' > and at n = nc, the charge susceptibil- 
ity diverges at a finite temperature which decreases by 
decreasing t' . As i' 0, ric moves continuously from 
values close to half filling to the filling that corresponds 
to the QCP. For t' = 0, the SC dome is centered at the 
QCP where the pseudogap temperature, T* , also van- 
ishes. This suggests that the incipient phase separation 
might play a role in the pairing mechanism. However, for 
t' > 0, while T* seems to vanish roughly at the center 
of the dome, phase separation happens at temperatures 
much higher than the SC temperature with classical criti- 
cal points that move to the lower-doping side of the dome. 
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an enhancement of pairing interactions^^—. Finally, it 
has been suggested that in the vicinity of the QCP, even 
a weak retarded attractive interaction may become far 
more effective at inducing pairing^. The relation be- 
tween the QCP and superconductivity will be explored 
in greater detail in future studies. 
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